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Abstract 

New CCD photometric observations of BX Dra were obtained for 26 nights from 
2009 April to 2010 June. The long-term photometric behaviors of the system are 
presented from detailed studies of the period and light variations, based on the his- 
torical data and our new observations. All available light curves display total eclipses 
at secondary minima and inverse O'Connell effects with Max I fainter than Max II, 
which are satisfactorily modeled by adding the slightly time-varying hot spot on the 
primary star. A total of 87 times of minimum light spanning over about 74 yrs, in- 
cluding our 22 timing measurements, were used for ephemeris computations. Detailed 
analysis of the O-C diagram showed that the orbital period has changed in combi- 
nations with an upward parabola and a sinusoidal variation. The continuous period 
increase with a rate of +5.65 xl0~ 7 d yr~ x is consistent with that calculated from 
the Wilson- Devinney synthesis code. It can be interpreted as a mass transfer from 
the secondary to the primary star at a rate of 2.74xl0 -7 M yr _1 , which is one of 
the largest rates for contact systems. The most likely explanation of the sinusoidal 
variation with a period of 30.2 yrs and a semi-amplitude of 0.0062 d is a light-travel- 
time effect due to the existence of a circumbinary object. We suggest that BX Dra is 
probably a triple system, consisting of a primary star with a spectral type of F0, its 
secondary component of spectral type Fl-2, and an unseen circumbinary object with 
a minimum mass of M3 = 0.23 M Q . 
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1. Introduction 

BX Dra (GSC 4192-0448, 2MASS J16061736+6245460, HIP 78891) was discovered to 
be a short-period variable by Strohmeier (1958) and classified as an RR Lyr-type star with 



a period of 0.561192 days by Strohmeier et al. (1965). Some doubt about this classification 
and hints about a binary nature were presented by Smith (1990), who proposed this star to 
be an ellipsoidal-type variable. Agerer & Dahm (1995) reclassified the variable as a (3 Lyr- 
type eclipsing binary from their photographic and CCD photometry without any filter. They 
also determined new linear ephemeris from the CCD measurements and suggested that all 
available timings could be represented by a quadratic ephemeris, in which the parabolic term 
(+5.56 xlO^ 10 d) indicates a continuous period increase. Pych et al. (2004) obtained double- 
line radial-velocity curves of BX Dra with semi-amplitudes of K\ = 80.0 km s _1 and Ki = 
276.4 km s _1 and found that this system is an A-subtype contact binary with a spectral type 
of F0IV-V. 

Recently, Sanchez-Bajo et al. (2007, hearafter SGG), Kim et al. (2009, hearafter KIM), 
and Zola et al. (2010, hearafter ZOLA) made CCD light curves in the BVI, BV, and BVRI 
bandpasses, respectively. Assuming both components to have a convective envelope and con- 
sidering a cool spot on the primary component, SGG analyzed their light curves and found a 
third light contributing 2.4% light in the B bandpass, 2.1% in V, and 0.9% in I. In addition, 
they computed the absolute dimensions of BX Dra from their photometric parameters and the 
spectroscopic orbit of Pych et al. (2004) and suggested the third light may come from a field 
object which is not bound to the eclipsing pair. On the other hand, ZOLA presented results of 
the modelling of their multicolor light curves without considering either a third light source or 
a spot. 

Although eclipsing minimum epochs have been reported assiduously by numerous work- 
ers, the period variation of BX Dra has not yet been studied in detail. In this article, we 
present and discuss the long-term photometric behaviors of the binary system together with 
the first detailed analysis of the O-C diagram, based on our new CCD observations as well as 
the historical data. 

2. CCD Photometric Observations 

New CCD photometric observations of BX Dra were obtained during two observing 
seasons in 2009 April and between 2010 April and June, using BVR filters attached to the 61- 
cm reflector at Sobaeksan Optical Astronomy Observatory (SOAO) in Korea. The observations 
of the 2009 season (SOAO09) were made on eight nights using a SITe 2K CCD camera, which 
has 2048x2048 pixels and an image field-of-view (FOV) of about 20'.5x20'.5 at the f/13.5 
Cassegrain focus of the telescope. The observations of the 2010 season (SOAO10) were made 
on 18 nights using an FLI IMG4301E 2K CCD camera. The new CCD chip has 2084x2084 
pixels and an FOV of about 20'.9x20'.9. The instruments and reduction methods were the 
same as those described by Lee et al. (2007, 2011). GSC 4192-0521 (V = 11.0) and GSC 
4192-0617 (V = 11.3) in the same observing field and with a color index similar to BX Dra 
were selected as comparison and check stars, respectively. We detected no variations in the 
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brightness difference between these two stars during our observing runs. An observed image in 
SOAO09 is given as Figure 1, in which the comparison and check stars are marked C and K, 
respectively. The la values of the dispersion of the (K— C) differences are about ±0.009 mag 
in all bandpasses. 

From the SOAO observations, we obtained 1620, 1622, and 1617 individual points in the 
B, V, and R bandpasses, respectively, and the sample is listed in Table 1. The light curves 
of BX Dra are plotted in Figure 2 as differential magnitude versus orbital phases, which were 
computed according to the ephemeris for our hot-spot model described in section 3. The open 
circles and plus symbols are the individual measures of the 2009 and 2010 seasons, respectively. 

3. Light-Curve Analysis And Spot Model 

As do historical light curves, our light curves of BX Dra show a flat bottom at the 
secondary minimum, which means that this system belongs to the A-subtype of W UMa-type 
stars. In addition, the new data show the inverse O'Connell effect with Max I (at phase 
0.25) fainter than Max II (at phase 0.75) by about 0.014 and 0.008 mag for the B and V, 
respectively, while the light levels in R are equal at the quadratures. These features usually 
indicate wavelength-dependent spot activity on the component stars. In order to obtain a 
unique set of photometric solutions, we analyzed all available light curves (SGG, ZOLA, KIM, 
SOAO) using contact mode 3 of the 2003 version of the Wilson-Devinney synthesis code (Wilson 
& Devinney 1971, hereafter W-D). For this purpose, we normalized the level at phase 0.75 and 
used a weighting scheme identical to that for the eclipsing binary GW Gem (Lee et al. 2009a). 
Table 2 lists the light-curve sets for BX Dra analyzed in this article and the standard deviations 
(a) of a single observation. 

The binary parameters of BX Dra were initialized in a manner similar to that for the 
contact systems BX Peg (Lee et al. 2004) and A A UMa (Lee et al. 2011). The effective tem- 
perature of the more massive primary star was fixed at T\ = 6,980 K, according to the spectral 
type F0 classified by Pych et al. (2004). Linear bolometric (X, Y) and monochromatic (x, y) 
limb-darkening coefficients were interpolated from the values of van Hamme (1993) in concert 
with the model atmosphere option. The initial value for the mass ratio (q = m 2 /m 1 ) was taken 
from Pych et al. (2004). Throughout the analyses, synchronous rotation was assumed and 
a third light source (£3) was considered. In addition, as the atmosphere of both components 
should lie close to the boundary between the radiative (hereafter RE) and convective (hereafter 
CE) envelopes, the gravity-darkening exponents (g) and the bolometric albedos (A) were inves- 
tigated at standard values of (A, g) = (1.0, 1.0) and (0.5, 0.32) for the two cases, respectively. 
Because the CE model gives a better fit than does the RE model, the common envelope was 
treated as a convective atmosphere in all subsequent syntheses. In Tables 3 and 4, parentheses 
represent adjusted parameters and the subscripts 1 and 2 refer to the primary and secondary 
stars being eclipsed at Min I and Min II, respectively. 
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First of all, we analyzed simultaneously all the light curves of BX Dra, permitting no 
spots. The unspotted solution is listed in the second column of Table 3 and the V residuals 
from the analysis are plotted in the left panels of Figure 3. Similar patterns exist for the 
other bandpasses. As shown in these panels, the model light curves do not fit the observed 
ones at all well. In contact binaries, the discrepancies may be caused by local photospheric 
inhomogeneities such as a cool spot on a magnetically active component and/or a hot spot due 
to impact from a mass transfer between the components (see, e.g., Lee et al. 2009b, 2010). 
Thus, spot models were added to fit the light variations. Because the seasonal variations of 
the light residuals are not large, we reanalyzed the datasets considering each of a hot and cool 
spot on each component. The results are given in columns (3)-(6) of Table 3 together with the 
spot parameters. The residuals from the hot-spot model (Hot 1) on the primary component 
are plotted in the right panels of Figure 3. From these results, we can see that the Hot 1 
model gives slightly smaller values for the sum of the residuals squared (T>W(0 — C) 2 ) than 
do the other models and that a large hot spot on the primary star has been sufficient for the 
light-curve representations of BX Dra. Nonetheless, it is difficult to distinguish between the 
spot models from only the light-curve analysis because the differences among them are small. 

Finally, to understand the long-term spot behaviors in detail, we solved five historical 
data sets separately using the hot-spot model on the primary star. For this procedure, we 
adjusted only the orbital ephemeris, spot, and luminosity parameters among the Hot 1 model 
parameters. The final results are given in Table 4 and the normalized V observations with the 
model light curves are plotted in Figure 4. The light ratios have not changed as a result of the 
spot modeling and the seasonal light curves could be represented by a slightly variable hot spot. 
Because the colatitudes and longitudes of the spot have been almost constant with time and 
the components of BX Dra should have shallow convective shells and at most weak magnetic 
activity as surmised from their temperatures, it is possible to regard the main cause of the 
activity as a sporadically-variable mass transfer from the secondary to the primary component. 
However, the streaming gas might usually lead to a hot spot around the primary equator, rather 
than the modeled colatitude very different from 90°. In reality, the intrinsic light variations of 
BX Dra could be caused by the simultaneous existence of a hot spot due to mass transfer and 
a magnetic spot on a component, but separate trials for such spot configurations did not give 
a better fit than the single hot-spot model. In all procedures, we looked for a possible third 
light source (£ 3 ) as suggested by SGG but found that the parameter remains indistinguishable 
from zero within its error. 

In order to measure the physical properties of the binary system, we analyzed the radial- 
velocity curves of Pych et al. (2004) with our light-curve parameters for the Hot 1 model. From 
the photometric and spectroscopic results, the absolute dimensions of BX Dra were determined 
and listed in Table 5, where the radii (R) are the mean- volume radii evaluated using tabulations 
of Mochnacki (1984). The luminosity (L) and bolometric magnitudes (Mbd) were computed 
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by adopting T e fj = 5,780 K and Mboi© = +4.73 for solar values. To estimate the uncertainty 
in the luminosity, it was assumed that the temperature of each component have an error of 
200 K in accordance with the unreliability in the spectral classification. For the absolute visual 
magnitudes (My), we used the bolometric corrections (BCs) from the scaling between logT 
and BC given by Torres (2010). The absolute parameters presented in this paper are consistent 
with those of SGG within the uncertainties. 

4. Orbital Period Study 

From the SOAO observations, 12 weighted times of minimum light and their errors were 
determined using the method of Kwee & van Woerden (1956). In addition to these, eight 
and two eclipse timings were newly determined from the individual measurements of KIM and 
ZOLA, respectively. Including our measurements, a total of 87 timings (36 photographic, 9 
photoelectric, and 42 CCD) have been collected from the database of Kreiner et al. (2001) and 
from more recent literature. All available photoelectric and CCD timings are listed in Table 6. 
Because many timings of the system have been published with no errors, the following standard 
deviations were assigned to timing residuals based on the observational method: ±0.0316 d for 
photographic plate, ±0.0019 d for photoelectric and ±0.0012 d for CCD minima. Relative 
weights for the period analysis of BX Dra were then scaled from the inverse squares of these 
values (Lee et al. 2007). 

As mentioned in the Introduction, SGG reported that the period change of BX Dra 
could be represented by a parabolic variation. After testing several other forms including this 
possibility, we found that the eclipse timings display a sinusoidal variation superposed on an 
upward parabola, rather than varying in a monotonic pattern. Using the PERIOD04 program 
(Lenz Sz Breger 2005), we looked to see if the residuals from the quadratic fit represent real 
and periodic variations. As shown in the small box of the top panel of Figure 5, a frequency 
of / = 0.0000636 cycle d" 1 was detected corresponding to about 25 yr. The periodic variation 
suggests a light-travel time (LTT) effect driven by the existence of a third component orbiting 
the eclipsing binary. Thus, the complete timing data were fitted to the following quadratic plus 
LTT ephemeris: 

C = T + PE + AE 2 + r 3 , (1) 

where T3 is the LTT due to a circumbinary object in the system (Irwin 1952, 1959) and includes 
five parameters (a^sin^, e, u, n and T). The Levenberg-Marquardt algorithm (Press et 
al. 1992) was applied to solve for the eight parameters of the ephemeris and the results are 
summarized in Table 7, together with related quantities. Our absolute dimensions given in 
Table 5 have been used for these and subsequent calculations. 

The O-C residuals calculated from the linear terms in equation (1) are plotted in the top 
panel of Figure 5, where the continuous curve represents the quadratic term of this ephemeris. 
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The middle panel displays the LTT orbit, and the bottom panel shows the residuals from the 
complete ephemeris. These appear as 0-Cf u n in the fourth column of Table 6. As displayed 
in the figure, the quadratic plus LTT ephemeris currently provides a good representation of all 
modern times of minimum light. The LTT orbit has a period of P3 = 30.2 yr, a semi-amplitude 
of K = 0.0062 d, a projected orbital semi-major axis of ai2sinz3 = 1.08 AU, and an eccentricity 
of e = 0.35. The mass function of the circumbinary object becomes f^(M^) = 0.00138 M Q 
and its minimum mass is 0.23 M Q . Because only about 62 % of the 30-year period has been 
covered by the photoelectric and CCD data, future accurate timings are required to identify 
and understand the LTT effect. 

As in the case of the contact binary AR Boo (Lee et al. 2009b) with a convective 
envelope, the quasi-sinusoidal period variation with small amplitudes may be produced by 
asymmetrical eclipse minima due to spot activity (Kalimeris et al. 2002) and/or the method 
of measuring the timings (Maceroni & van't Veer 1994). The light-curve synthesis method 
gives more and better information with respect to the other methods, which do not consider 
spot activity and are based only on observations during minima (Lee et al. 2009b). Because 
five datasets of BX Dra were modeled for Hot 1 spot parameters, we calculated a minimum 
epoch for each eclipse in these datasets with the W-D code by means of adjusting only the 
ephemeris epoch (T ). The results are listed in Table 8 together with the previously tabulated 
timings for comparison and are illustrated with the 'x' symbols in Figure 5. The differences 
between the published minima and those obtained by the synthesis method are significantly 
smaller than the observed amplitude (0.012 d) of the LTT variation. As shown in Figure 5, 
the light-curve timings agree with our analysis of the O—C diagram and the periodic variation 
cannot result from the starspot activity. Further, as shown in the fourth column of Table 8, 
there are systematic runs of differences between them, which are negative for Min I and positive 
for Min II. These differences are caused by the hot spot on the primary star presented to the 
observer. 

The quadratic term (A) of equation (1) signifies a continuous period increase with 
a rate of + (5.65±0.07)xl0 -7 d yr" 1 , corresponding to a fractional period change of 
+ (1.55±0.02) xlO" 9 . This is very close to the value of +1.29xl0~ 9 calculated from the W-D 
code, independently of the eclipse timings. The most common explanation of the period in- 
crease is a mass transfer from the less massive secondary star to the primary component. Under 
the assumption of conservative mass transfer, the transfer rate is calculated to be 2.74xl0~ 7 
M yr" 1 . This value is larger than those recently derived for other overcontact systems [e.g., 
2.0 xlO" 8 for BV Dra (Yang et al. 2009), 1.5 xlO" 7 for AR Boo (Lee et al. 2009b), 7.4 xlO" 8 
for V1191 Cyg (Zhu et al. 2011), and 6.6xl0" 8 for AA UMa (Lee et al. 2011), all in M yr" 1 
units] . 
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5. Summary And Discussion 

In this article, we have presented the long-term photometric behaviors of BX Dra from 
the detailed analyses of the light curves and the O-C diagram based on the historical and new 
observations. The results from these analyses can be summarized as follows: 

1. Historical light curves of BX Dra, as well as our own, display total eclipses at secondary 
minima and inverse O'Connell effects with Max I fainter than Max II. The asymmetric 
light curves can satisfactorily be explained by the spot model. The slightly variable hot 
spot on the primary star permits good light-curve representations for the system. 

2. The orbital period of BX Dra has varied in a combination with an upward parabola and 
a sinusoidal variation with a period of 30.2 yrs and a semi-amplitude of 0.0062 d, rather 
than in a monotonic fashion. The increasing rate of the secular period is calculated to be 
+5.65xl0~ 7 d yr -1 , which may be caused by the mass transfer from the secondary to the 
primary component in the system with a rate of about 2.74 xlO -7 M Q yr _1 . 

3. The sinusoidal variation can be interpreted as the LTT effect due to the existence of a 
circumbinary object. If the third companion is on the main sequence and its orbit is 
coplanar with the eclipsing pair (i 3 = 81°. 8), the mass of the circumbinary object is M 3 
= 0.23 M and its radius and temperature are calculated to be R3 = 0.24 R and T 3 
= 3022 K, respectively, following the empirical relations from Southworth (2009). These 
correspond to a spectral type of about M6 V and a bolometric luminosity of L 3 = 0.004 L Q 
and contribute about 0.03% to the total light of the triple system. So, it would be difficult 
to detect such a faint companion from analyses of spotted light curves and spectroscopic 
observations. The absence of the evidence does not rule out the presence of the hypothetic 
third component. 

4. The results presented in this paper indicate that BX Dra is an A-subtype overcontact 
binary with an unseen circumbinary object; the more massive primary star has a spectral 
type of F0 and the secondary component a spectral type of Fl-2. We think that the hot 
spot on the primary star could be produced as a result of the impact of the gas stream from 
the less massive secondary star. The higher mass transfer rate than the other overcontact 
systems and no variation of spot parameters for about four years support this hot spot 
model. Moreover, because both components should not have deep convective envelopes as 
surmised from their temperatures, it is not reasonable to imagine that magnetic spots on 
the components of the system is the dominant cause of the light variation. 

Although all historical timings of BX Dra are in good agreement with the calculated LTT 
effect as seen in Figure 5, because of the absence of any independent third body detection, we 
consider the possibility that a magnetic activity cycle may be the main cause of the sinusoidal 
period modulation (Applegate 1992, Lanza et al. 1998). With the period (P 3 ) and amplitudes 
(K) listed in Table 7, the model parameters were calculated for each component from the 
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Applegate formulae. The parameters are listed in Table 9, where the rms luminosity changes 
(Am rms ) converted to magnitude scale were obtained with equation (4) in the paper of Kim 
et al. (1997). The variations of the gravitational quadrupole moment (AQ) are two orders 
of magnitude smaller than typical values of 10 51 ~ 10 52 for close binaries (Lanza & Rodono 
1999). In reality, the magnetic mechanism is adequate for systems with a spectral type later 
than about F5 (Hall 1989), unlike the case of BX Dra. Moreover, it is difficult for this model 
to produce perfectly smooth and tilted periodic components in the O-C variation. At present, 
because there is no other alternative but the LTT effect, the sinusoidal variation most likely 
arises from an unseen third companion gravitationally bound to the eclipsing binary. The 
circumbinary object in BX Dra may have played an important role in the formation of the 
eclipsing pair, which would finally evolve into single star by angular momentum loss through 
magnetic braking. Future high-precision timing measurements will be crucial unveiling the 
orbital period change of this system. 

We would like to thank the staff of the Sobaeksan Optical Astronomy Observatory for 
assistance during our observations. We also thanks F. Agerer, F. Sanchez-Bajo, and S. Zola 
for sending us their published data on BX Dra. This research has made use of the Simbad 
database maintained at CDS, Strasbourg, France. This work was supported by the KASI 
(Korea Astronomy and Space Science Institute) grant. 
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BX Dra 



•t - 



Fig. 1. A sample of observed CCD images (20'.5x20'.5) of BX Dra from the 2009 season. Monitoring 
numerous frames led us to choose star C as a comparison and K as a check. North is up and east is to the 
left. 
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Fig. 2. BVR light curves of BX Dra observed in 2009 and 2010. Because of the high density of the points, 
many of the 2009 measures cannot be seen individually. 
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Unspotted Model 



Hot 1 Model 




Fig. 3. Left and right panels show the light residuals from the solutions without a spot and with a hot 
spot (Hot 1) listed in Table 3, respectively. The data refer to the V bandpass. 
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Fig. 4. Normalized V observations with the fitted model light curves. The continuous curves represent 
the solutions obtained with the hot-spot model parameters listed in Table 4. 
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Fig. 5. O-C diagram of BX Dra. In the top panel, constructed with the linear terms of the quadratic 
plus LTT ephcmeris, the continuous curve represents the quadratic term of the equation. The inset box 
is drawn with the result from the PERIOD04 formalism for the residuals, where a dominant frequency of 
/ = 0.0000636 cycle d _1 is found with a semi-amplitude of 0.0054 mag. The middle panel displays the 
LTT orbit and the bottom panel shows the CC and PE residuals from the complete ephcmeris. The 'x' 
symbols refer to new minimum timings obtained from the W-D code. 
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Table 1. CCD photometric observations of BX Dra in 2009 and 2010.* 
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-1 


.110 


2,454,927. 


.18680 


-0 


.764 


2,454,927.18545 


-0 


.573 










2010 












2,455,303.26231 


-1 


.065 


2,455,303. 


26331 


-0 


.724 


2,455,303.26406 


-0 


.521 


2,455,303.26501 


-1 


.077 


2,455,303. 


26602 


-0 


.742 


2,455,303.26677 


-0 


.529 


2,455,303.26770 


-1 


.084 


2,455,303. 


.26872 


-0 


.749 


2,455,303.26946 


-0 


.547 


2,455,303.27049 


-1 


.105 


2,455,303. 


27163 


-0 


.774 


2,455,303.27245 


-0 


,558 


2,455,303.27351 


-1 


.125 


2,455,303. 


27465 


-0 


.792 


2,455,303.27548 


-0 


.588 



* A sample is shown here: the full version is available in its entirety 
in a machine-readable form in the PASJ online edition. 
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Table 2. Light-curve sets for BX Dra. 



Reference 


Season 


Data Type 




SGG 


2006 


B 


0.0078 






V 


0.0062 






I 


0.0058 


ZOLA 


2006 


B 


0.0115 






V 


0.0121 






R 


0.0120 






I 


0.0113 


KIM 


2008 


B 


0.0209 






V 


0.0119 




900Q 


n 
u 


U.UU ( 






V 


0.0087 






R 


0.0093 




2010 


B 


0.0098 






V 


0.0070 






R 


0.0062 



a In units of total light at phase 0.75. 
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Table 3. Binary parameters obtained by fitting simultaneously all light curves. 



Parameter 


Without Spot 




Spot Model a 








Cool 1 


Cool 2 


Hot 1 


Hot 2 


To (HJD) b 


810.58987(12) 


810.59009(13) 


810.58995(13) 


810.58976(12) 


810.58970(12) 


P(d) 


0.57902455(4) 


0.57902476(4) 


0.57902477(4) 


0.57902479(4) 


0.57902472(4) 


dP/dt (xlO" 9 ) 


1.386(12) 


1.303(11) 


1.302(11) 


1.288(11) 


1.317(11) 


T X (K) 


6980 


6980 


6980 


6980 


6980 


T2 (K) 


6979(2) 


6707(2) 


6995(2) 


6758(4) 


6805(3) 


i (deg) 


80.63(6) 


80.66(5) 


82.42(8) 


81.80(6) 


82.18(8) 


q ( = mj/mi) 


0.2884(5) 


0.2874(5) 


0.2871(4) 


0.2884(5) 


0.2882(5) 


fii = 2 


2.3475(16) 


2.3462(15) 


2.3404(14) 


2.3308(14) 


2.3293(16) 


Li/(Li+i 2 )s 


0.7432(5) 


0.7828(6) 


0.7407(4) 


0.7729(4) 


0.7660(5) 


il/(il + L 2 )y 


0.7435(5) 


0.7744(5) 


0.7416(4) 


0.7660(4) 


0.7606(4) 


L 1 /(L 1 +L 2 )r 


0.7437(4) 


0.7685(5) 


0.7421(4) 


0.7612(4) 


0.7568(4) 


L 1 /(L 1 +L 2 ) I 


0.7437(4) 


0.7639(4) 


0.7424(4) 


0.7575(3) 


0.7539(4) 


n (pole) 


0.4791(4) 


0.4792(4) 


0.4804(3) 


0.4828(3) 


0.4832(4) 


ri (side) 


0.5218(5) 


0.5219(5) 


0.5237(5) 


0.5272(5) 


0.5277(5) 


ri (back) 


0.5542(7) 


0.5543(7) 


0.5565(6) 


0.5614(7) 


0.5620(7) 


r 2 (pole) 


0.2803(6) 


0.2799(6) 


0.2811(5) 


0.2846(6) 


0.2849(6) 


r 2 (side) 


0.2953(8) 


0.2949(7) 


0.2964(7) 


0.3007(8) 


0.3010(8) 


r 2 (back) 


0.3513(18) 


0.3506(17) 


0.3542(16) 


0.3638(20) 


0.3648(21) 


Spot parameters: 












Colatitudei (deg) 




22.3 




8.5 




Longitudei (deg) 




188.8 




11.5 




Radiusi (deg) 




35.7 




36.5 




Tspotl/Tiocall 




0.758 




1.361 




Colatitude 2 (deg) 






24.1 




33.3 


Longitude 2 (deg) 






9.6 




186.0 


Radius 2 (deg) 






50.4 




33.2 


T S pot2/T'iocal2 






0.723 




1.259 


T,W{0-C) 2 


0.0151 


0.0129 


0.0126 


0.0121 


0.0126 


a Cool 1: a cool spot on the primary; Cool 2: a cool spot on the secondary; Hot 1: a hot spot 


on the primary; Hot 2: a hot spot 


on the secondary. 









b HJD 2,449,000 is suppressed. 
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Table 4. Year-to-year variations of the spot and luminosity parameters. 



Parameter 


SGG 


ZOLA 


KIM 


SOAO09 


SOAO10 


T (HJD)° 


3767.66365(7) 


3905.47131(5) 


4581.20182(19) 


4932.09274(8) 


5314.25365(7) 


P(d) 


0.5790304(16) 


0.5790455(73) 


0.5790273(20) 


0.5790446(45) 


0.5790425(13) 


Colatitudei (cleg) 


8.5 


8.5 


8.8 


8.5 


8.3 


Longitudei (deg) 


10.9 


11.5 


11.5 


11.5 


11.5 


Radiusx (deg) 


36.1 


39.1 


38.5 


39.6 


36.5 


T IT 

spotX/ loceilX 


1.365 


1.369 


1.386 


1.340 


1.370 


L 1 /(L 1 +L 2 ) B 


0.7729(4) 


0.7729(3) 


0.7729(10) 


0.7729(4) 


0.7729(2) 


L 1 /(L 1 +L 2 ) V 


0.7660(3) 


0.7660(3) 


0.7660(5) 


0.7660(3) 


0.7660(2) 


L 1 /{L 1 +L 2 )r 




0.7612(3) 




0.7612(3) 


0.7612(2) 


L 1 /{L 1 +L 2 )i 


0.7575(2) 


0.7575(3) 








XW(0-C) 2 


0.0092 


0.0119 


0.0113 


0.0121 


0.0116 



a HJD 2,450,000 is suppressed. 
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Table 5. Absolute parameters for BX Dra. 



Parameter 


Primary 


Secondary 


a (R Q ) 


4.058(87) 


V a (km s- 1 ) 


-25.6(2.2) 


M (Mq) 


2.08(10) 


0.60(4) 


R (R Q ) 


2.13(5) 


1.28(3) 


log g (cgs) 


4.10(3) 


4.00(3) 


p (g. cm -3 ) 


0.30(2) 


0.41(4) 


T(K) 


6980(200) 


6758(200) 


L (i ) 


9.66(1.18) 


3.05(38) 


Mboi (mag) 


2.27(13) 


3.52(14) 


BC (mag) 


0.03 


0.02 


M v (mag) 


2.24(13) 


3.50(14) 
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Table 6. Photoelectric and CCD timings of minimum light for BX Dra. 



HJD 

(2,400,000+) 


Error 


Epoch 


o-cw 


Min 


References 


48,528.63 




-2214 





0.0021 


I 


Agerer & Dahm (1995) 


49,810.5926 










0.0008 


I 


Agerer & Dahm (1995) 


49,811.4614 




1 


5 


0.0010 


11 


Agerer & Dahm (1995) 


49,812.3275 




3 





-0.0014 


I 


Agerer & Dahm (1995) 


49,840.4122 




51 


5 


0.0005 


11 


Agerer & Dahm (1995) 


49,866.4682 




96 


5 


0.0002 


II 


Agerer & Dahm (1995) 


49,888.4720 




134 


5 


0.0010 


II 


Agerer & Dahm (1995) 


50,147.5838 


±0.0003 


582 





-0.0019 


I 


Agerer & Hiibscher (1996) 


50,547.4042 


±0.0004 


1272 


5 


0.0003 


11 


Agerer & Hiibscher (1999) 


50,904.3743 


±0.0003 


1889 





0.0004 


I 


Agerer & Hiibscher (1999) 


50,945.4868 


±0.0005 


1960 





0.0020 


I 


Agerer & Hiibscher (1999) 


51,256.4227 


±0.0004 


2497 





0.0007 


I 


Agerer & Hiibscher (2000) 


51,270.6090 


±0.0011 


2521 


5 


0.0009 


II 


Agerer & Hiibscher (2000) 


52,362.360 


±0.003 


4407 





-0.0025 


I 


Diethelm (2002) 


52,401.4508 


±0.0008 


4474 


5 


0.0039 


11 


Agerer & Hiibscher (2003) 


53,409.5346 


±0.0007 


6215 


5 


0.0005 


11 


Hiibscher ct al. (2005) 


53,758.6891 


±0.0002 


6818 


5 


0.0006 


II 


Sanchcz-Bajo ct al. (2007) 


53,767.6644 


±0.0003 


6834 





0.0010 


I 


Sanchcz-Bajo ct al. (2007) 


53,772.5851 


±0.0001 


6842 


5 


-0.0002 


II 


Sanchez-Bajo et al. (2007) 


53,774.6127 


±0.0001 


6846 





0.0009 


I 


Sanchcz-Bajo ct al. (2007) 


53,794.5882 


±0.0001 


6880 


5 


-0.0001 


11 


Sanchcz-Bajo ct al. (2007) 


53,800.6687 


±0.0001 


6891 





0.0006 


I 


Sanchcz-Bajo ct al. (2007) 


53,803.5635 


±0.0001 


6896 





0.0002 


I 


Sanchcz-Bajo ct al. (2007) 


53,829.9089 


±0.0002 


6941 


5 


-0.0002 


II 


Nelson (2007) 


53,846.4111 


±0.0004 


6970 





-0.0003 


I 


Diethelm (2006) 


53,905.4719 


±0.0002 


7072 





-0.0006 


I 


Zola ct al. (2010) 


53,907.4969 


±0.0002 


7075 


5 


-0.0021 


II 


Zola et al. (2010) 


54,527.9296 


±0.0001 


8147 





-0.0001 


1 


Nelson (2009) 


54,575.1205 


±0.0002 


8228 


5 


-0.0002 


II 


Kim ct al. (2009) 


54,581.2027 


±0.0003 


8239 





0.0021 


I 


Kim et al. (2009) 


54,584.0982 


±0.0006 


8244 





0.0025 


I 


Kim et al. (2009) 


54,586.1232 


±0.0007 


8247 


5 


0.0009 


II 


Kim et al. (2009) 


54,588.1501 


±0.0005 


8251 





0.0012 


I 


Kim ct al. (2009) 


54,597.4133 


±0.0001 


8267 





-0.0001 


I 


Hiibscher (2009) 


54,603.2050 


±0.0006 


8277 





0.0013 


I 


Kim ct al. (2009) 


54,685.1390 


±0.0003 


8418 


5 


0.0024 


11 


Kim ct al. (2009) 


54,686.0070 


±0.0005 


8420 





0.0018 


I 


Kim ct al. (2009) 


54,931.2241 


±0.0003 


8843 


5 


-0.0007 


II 


This paper (SOAO) 


54,932.0935 


±0.0003 


8845 





0.0001 


I 


This paper (SOAO) 


54,934.1202 


±0.0003 


8848 


5 


0.0002 


II 


This paper (SOAO) 


54,941.3595 


±0.0003 


8861 





0.0016 


I 


Brat et al. (2009) 


54,951.2017 


±0.0002 


8878 





0.0003 


I 


This paper (SOAO) 


54,955.8323 


±0.0012 


8886 





-0.0014 


I 


Diethelm (2009) 


55,314.2540 


±0.0004 


9505 





-0.0005 


I 


This paper (SOAO) 


55,327.2826 


±0.0002 


9527 


5 


-0.0000 


11 


This paper (SOAO) 


55,332.2046 


±0.0001 


9536 





0.0001 


I 


This paper (SOAO) 


55,337.1261 


±0.0002 


9544 


5 


-0.0002 


11 


This paper (SOAO) 


55,350.1545 


±0.0002 


9567 





0.0000 


I 


This paper (SOAO) 


55,352.1809 


±0.0002 


9570 


5 


-0.0002 


II 


This paper (SOAO) 


55,353.0485 


±0.0002 


9572 





-0.0012 


I 


This paper (SOAO) 


55,354.2074 


±0.0003 


9574 





-0.0002 


I 


This paper (SOAO) 
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Table 7. Parameters for the quadratic plus LTT ephemeris of BX Dra. 



Parameter 


Values 


Unit 


To 


2,449,810.58844(36) 


HJD 


P 


0.579024741(46) 


(1 


A 


4.476(56) xl0~ 10 


d 


ai 2 sirn 3 


1.08(13) 


AU 


w 


61.8(3.8) 


deg 


e 


0.350(94) 




n 


0.0326(13) 


deg d _1 


T 


2,450,417(140) 


HJD 


P 3 


30.2(1.2) 


yr 


K 


0.00615(74) 


d 


/(Ms) 


0.00138(17) 


M Q 


M3 sini3 


0.23 


M q 


dP/dt 


5.647(70) xl0~ 7 


d yr -1 


dM/dt 


2.741 xl0~ 7 


M Q yr' 1 
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Table 8. Minimum timings determined by the W-D code from individual eclipses of BX Dra. 



Observed a ' b 


W-D b 


Error 


Diffcrcnce d 


Filter 


Min 


References 




3,758.6891 


3,758.68907 


±0.00023 


0.00000 


BV 


II 


Sanchcz-Bajo et al. 


(2007) 


3,767.6644 


3,767.66375 


±0.00015 


+0.00067 


BV 


I 


Sanchcz-Bajo et al. 


(2007) 


3,772.5851 


3,772.58546 


±0.00022 


-0.00041 


BV 


II 


Sanchcz-Bajo et al. 


(2007) 


3,774.6127 


3,774.61214 


±0.00016 


+0.00056 


BV 


I 


Sanchcz-Bajo et al. 


(2007) 


3,794.5882 


3,794.58853 


±0.00014 


-0.00032 


I 


II 


Sanchcz-Bajo et al. 


(2007) 


3,800.6687 


3,800.66844 


±0.00013 


+0.00029 


I 


I 


Sanchcz-Bajo et al. 


(2007) 


3,803.5635 


3,803.56325 


±0.00007 


+0.00023 


I 


I 


Sanchcz-Bajo et al. 


(2007) 


3,905.4719 


3,905.47111 


±0.00011 


+0.00075 


BVRI 


I 


Zola et al. (2010) 




3,907.4969 


3,907.49692 


±0.00014 


0.00000 


BVRI 


II 


Zola et al. (2010) 




4,575.1205 


4,575.12128 


±0.00047 


-0.00078 


BV 


II 


Kim et al. (2009) 




4,581.2027 


4,581.20221 


±0.00042 


+0.00049 


BV 


I 


Kim et al. (2009) 




4,584.0982 


4,584.09820 


±0.00047 


0.00000 


BV 


I 


Kim et al. (2009) 




4,586.1232 


4,586.12320 


±0.00076 


0.00000 


BV 


II 


Kim et al. (2009) 




4,588.1501 


4,588.14872 


±0.00053 


+0.00138 


BV 


I 


Kim et al. (2009) 




4,603.2050 


4,603.20341 


±0.00055 


+0.00159 


BV 


I 


Kim et al. (2009) 




4,685.1390 


4,685.13900 


±0.00047 


0.00000 


BV 


II 


Kim et al. (2009) 




4,686.0070 


4,686.00700 


±0.00070 


0.00000 


BV 


I 


Kim et al. (2009) 




4,931.2241 


4,931.22464 


±0.00028 


-0.00051 


BVR 


II 


This article 




4,932.0935 


4,932.09287 


±0.00018 


+0.00065 


BVR 


I 


This article 




4,934.1202 


4,934.12040 


±0.00014 


-0.00019 


BVR 


II 


This article 




4,951.2017 


4,951.20127 


±0.00018 


+0.00047 


BVR 


I 


This article 




5,314.2540 


5,314.25309 


±0.00017 


+0.00089 


BVR 


I 


This article 




5,327.2826 


5,327.28296 


±0.00018 


-0.00032 


BVR 


II 


This article 




5,332.2046 


5,332.20401 


±0.00009 


+0.00056 


BVR 


I 


This article 




5,337.1261 


5,337.12605 


±0.00018 


0.00000 


BVR 


II 


This article 




5,350.1545 


5,350.15382 


±0.00011 


+0.00063 


BVR 


I 


This article 




5,352.1809 


5,352.18145 


±0.00012 


-0.00058 


BVR 


II 


This article 




5,353.0485 


5,353.04813 


±0.00013 


+0.00033 


BVR 


I 


This article 




5,354.2074 


5,354.20687 


±0.00020 


+0.00056 


BVR 


I 


This article 





a cf. Table 6. 

b HJD 2,450,000 is suppressed. 

c Uncertainties yielded by the W-D code. 

d Differences between columns (1) and (2). 
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Table 9. Applcgatc parameters for possible magnetic activity of BX Dra. 



Parameter 


Primary 


Secondary 


Unit 


AP 


0.1753 


0.1753 


s 


AP/P 


3.50 x 10- 6 


3.50 x 10~ 6 




AQ 


1.28 x 10 50 


3.70 x 10 49 


g cm 2 


AJ 


2.60 x 10 47 


1.00 x 10 47 


? — i 
g cm s 


Is 


6.05 x 10 54 


6.30 x 10 53 


g cm 2 


An 


4.30 x 10~ 8 


1.58 x 10~ 7 


s' 1 


Ari/n 


3.42 x 10~ 4 


1.26 x 10" 3 




AE 


2.24 x 10 40 


3.15 x 10 40 


erg 




7.37 x 10 31 


1.04 x 10 32 


erg s~ 4 




0.0192 


0.0271 


L 




0.0020 


0.0089 


Ll,2 




±0.0016 


±0.0023 


mag 


B 


4.1 


5.4 


kG 
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